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Abstract: This paper presents a method of fusion of identification (attribute) information provided by two types 
of sensors: combined primary and secondary (IFF) surveillance radars and ESM (Electronic Support Measures). 
In the first section, the basic taxonomy of attribute identification is adopted in accordance with the standards 
of STANAG 1241 ed. 5 and STANAG 1241 ed. 6 (draft). These standards provide the following basic values of 
the attribute identifications: FRIEND, HOSTILE, NEUTRAL, UNKNOWN and additional values: ASSUMED 
FRIEND and SUSPECT. The basis of theoretical considerations is Dezert-Smarandache theory (DSmT) of 
inference. The paper presents and practically uses for combining identification information from different ESM 
sensors and radars six information fusion rules proposed by DSmT - the Proportional Conflict Redistribution 
rules (PCR1, PCR2, PCR3, PCR4, PCR5 and PCR6). In the paper, rules of determining attribute information by 
ESM sensor equipped with the data base of radar emitters are presented. It was proposed that each signal 
vector sent by the ESM sensor contained an extension specifying a randomized identification declaration 
(hypothesis) - basic belief assignment (BBA). The paper also presents a model for determining the basic belief 
assignment for a combined primary and secondary radar. Results of the PCR rules of sensor information 
combining for different scenarios of radio-electronic situation (deterministic and Monte Carlo) are presented 
in the final part of the paper. They confirm the legitimacy of the use of Dezert-Smarandache theory into 
information fusion for primary radars, secondary radars and ESM sensors. 


Keywords: information fusion; Dezert-Smarandache theory (DSmT) of inference; conflict 
redistribution rules; radar emitters recognition; electronic support measures (ESM); primary and 
secondary radars 


1. Introduction 


The paper is devoted to the fusion of identification information from ESM sensors and combined 
primary and secondary radar (IFF) using the rules of Dezert-Smarandach theory (DSmT) called 
proportional conflict redistribution rules. 

The first part of the paper presents the applied interpretation of attribute identification in 
accordance with the NATO STANAG 1241 standard. It should be noted that this is one of the possible 
interpretations of the adopted definitions. It leads to the Bayesian model of the basic belief 
assignment. 

The identification classification method depends on the organization that operates the ESM 
sensors. In the paper, one assumes that the sensor identification classification is consistent with the 
NATO STANAG 1241 standard [1,2]. In addition, one assumes that five identification classes are used 
- three primary and two secondary ones. Sensors can transmit identification information in the form 
of a hard decision, sometimes determined as non-randomized, or a soft decision, sometimes 
determined as a randomized decision. In the paper, one assumes that the sensors send identification 
information to the system in a randomized form, i.e. in the form of basic belief assignment on the set 
of identification classes. This assignment determines the sensor's belief that the detected emitter 
belongs to separate identification classes. 
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The next part of the paper presents the mathematical form of the DSmT conflict proportional 
redistribution rules PCR1, PCR2, PCR3, PCR4, PCR5 and PCR6 [3, 4] for two sensor inputs and PCR5 
and PCR6 for three sensor inputs, assuming the Bayesian model of the basic belief assignment of 
hypothesis. 

The next two sections show how to determine the basic belief assignment for combined primary 
and secondary (IFF) radar and ESM sensors. 

Combined primary and secondary (IFF) radars are the main source of identification information 
about air and maritime objects. A primary radar allows only to detect an object in a supervised area. 
The detection of the object is the precondition for sending a request to the object by the secondary 
radar (interrogator). Interpretation of the object response is dependent on the type of request. The so- 
called civilian modes allow only to determine whether the detected object replies to an interrogation 
or not. The paper presents a method for determining the basic belief assignment of airborne targets 
moving in observation space of combined primary and secondary (IFF) radars sensor. 

ESM (electronic intelligence - electronic support measures) electronic surveillance sensors 
consist of passive receivers and direction-finders, which allows them to capture emitter signals 
coming from certain directions. In this way, the electronic recognition system can receive, among 
others, information on radar emitters mounted on air or maritime platforms. Reports sent from the 
ESM sensors include, among others, the characteristics of the intercepted signal, the emitter's azimuth 
and the so-called identification information. 

The paper also assumes that sensors are equipped with specialized databases called the 
databases of emitter signal patterns, in which information about previously captured, processed, 
analyzed, recognized and described radar emitter signals is stored along with additional information 
about the type and mode of the emitter work, the platform on which these emitters can be installed, 
and the national or organizational affiliation of these platforms. The detected signals are the subject 
of an analysis procedure, which allows to determine the so-called distinctive features of the signal 
and then assigning this information to a specific electronic entity (already existing or created ad hoc) 
[5]. The basis for assigning distinctive information to an electronic entity is the azimuth angle of the 
incoming signal. 

In the case of a high density of targets, identification information may fluctuate due to incorrect 
assignment of signal information to the electronic entity [6]. The impact of this negative phenomenon 
can be significantly reduced by an efficient estimation of the emitter positions [7]. Assuming that 
sensors send all reports on the tracked electronic entities to the superior operation center in the 
electronic recognition system, such a center (in the paper called the information fusion center (IFC) 
can perform the fusion function of the identification information. The fusion of identification 
information ensures the greater stability of this information - resistance to accidental changes in 
sensor decisions. 

Each sensor report sent to the fusion information center contains a vector of belief mass for all 
attribute identification values. Results of the Proportional Conflict Redistribution sensor information 
combining rules for selected deterministic and Monte Carlo scenarios are presented in the final part 
of the paper. The identification information fusion can be realized based on three basic theories - 
Bayesian theory of inference, Dempster-Shafer theory - called the theory of evidence and Dezert- 
Smarandache theory. The methods of Dezert-Smarandache information fusion are used in this paper. 
In addition, their effectiveness is compared with the Dempster's rule of inference. 

At the end of the paper conclusions are given. They confirm the legitimacy of the use of Dezert- 
Smarandache theory into information fusion for primary radars, secondary radars and ESM sensors. 


2. Interpretation of attribute identification according to STANAG 1241 


The set of possible values of attribute identifications used by sensors can be adopted based on 
standardization documents of organizations that exploit these sensors [1,2,8,9,10]. 

This paper assumes a basic taxonomy of identification in accordance with the draft of STANAG 
1241 ed. 6 [2]. To other similar documents one may include the following standards: STANAG 4420 
and STANAG 1241 ed. 5, which provide the following basic values of the attribute identifications: 
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e FRIEND (F), 

e HOSTILE (H), 

e NEUTRAL (N), 

e UNKNOWN (UV). 

Each of these documents contain their own definitions of the declarations. 
The following definitions of these basic values of the attribute identification are used in the paper 

(in accordance with [2]): 

e FRIEND - an allied/coalition military track, object or entity; a track, object or entity, supporting 
friendly forces and belonging to an allied/coalition nation or a declared or recognized friendly 
faction or group, 

e HOSTILE - a track, object or entity whose characteristics, behavior or origin indicate that it 
belongs to opposing forces or poses a threat to friendly forces or their mission, 

e NEUTRAL - a military or civilian track, object or entity, neither belonging to allied/coalition 
military forces nor to opposing military forces, whose characteristics, behavior, origin or 
nationality indicates that it is neither supporting nor opposing friendly forces or their mission, 

e UNKNOWN -an evaluated track, object or entity, which does not meet the criteria for any other 
standard identity. 

These standards bring additional values of the attribute identification: 

e ASSUMED FRIEND, 

e SUSPECT. 

One should pay attention on these two recent identities contained in [1] as well as their 

definitions [2]: 

e ASSUMED FRIEND - a track, object or entity which is assumed to be friend or neutral because 
of its characteristics, behavior or origin, 

e SUSPECT - a track, object or entity whose characteristics, behavior or origin indicate that it 
potentially belongs to opposing forces or potentially poses a threat to friendly forces or their 
mission. 

The identification definitions in [1,2] can lead to different interpretations. This paper adopts the 

interpretation, the graphical form of which is shown in Figure 1. 


Figure 1. The interpretation of STANAG 1241 using the Venn diagram. 
3. Fusion of information from the ESM sensors and radars in the information fusion center (IFC) 


3.1. Diagram of the process of information fusion for two sensors in the information fusion center 


In this work, it is assumed that ESM sensors send messages asynchronously to the information 
fusion center. These reports contain sensor decisions regarding the identification of objects emitting 
detected signals. The set of possible identifications is following: 


0 = {6,,=1,...,6} (1) 
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wherein the following interpretation is used: 

@1 - FRIEND (F), 

62 - HOSTILE (H), 

@3 - NEUTRAL (N), 

64 - ASSUMED FRIEND(AF), 

@s - SUSPECT(S). 

06 - UNKNOWN (U), 

According to Figure 1, the hypotheses are mutually exclusive, i.e. 
6;, ifi = j, 


0,0 ={ D ifi xj. (2) 


Each sensor with the number i (i € N) sends its decisions as so-called soft decisions, i.e. as BBA 
measure vectors (BBA - basic belief assignment) 


m; = [m;(4,), ...,m;(46)]. (3) 
One should also introduce a vector of generalized BBA measures for the information fusion 
center 
mp = [mp( 44), ...,Mp( O6)]. (4) 


The paper adopts the Bayesian BBA model due to the fact that this model has been adopted as 
valid in the STANAG 4162 standard [9]. This means that equation (5) applies in addition to (1) and 


(2). 
Dik1 Mp(O;) = LP. m,( 6;) = 1. (5) 


The first case will be considered when two sensors send, asynchronously in one cycle, one report 
each containing decisions regarding the BBA related to the target. The IFC system, after receiving the 
report from the sensor, fuses the information contained in the two vectors: in the current generalized 
BBA vector m,; = [mp(6,),..,Mp(9.)] , and in the BBA vector m, from sensor 1 or in the BBA 
vector m, from sensor 2. 

The information fusion procedure performed in the IFC is carried out in accordance with the 
following formula: 


my; = Re(m,,m,) (i = 1 or 2), (6) 


wherein m;, is a vector of the generalized BBA measure determined by the R; rule based on 
the previous generalized BBA measure vector m; and the new BBA measure vector m; sent by the 
i-th sensor. The diagram of the identification information fusion from the ESM sensors is shown in 
Figure 2. 


Sensor 1 IFC 
Procedure of 


ESM sensor 
and radar 
report fusion 
for two bba 


vectors 


Figure 1. The diagram of the information fusion process in the information fusion center IFC for two 
sensors. Explanations: mi —- BBA measure vector of i-th sensor, mF — generalized BBA measure vector 
that is a part of the electronic entity record in IFC, EER - electronic entity record in IFC database. 
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The second case will be considered when two sensors send, asynchronously in one cycle, one 
report each containing decisions regarding the BBA related to the target. The IFC system waits for 
reports from both sensors in one cycle, using registers. Only when both registers are full, the IFC 
system performs a fusion of the information contained in three vectors: BBA vector m; = 
[mp(6,),-.,Mp(,)] , and BBA vector m, from sensor 1, and BBA vector m, from sensor 2. It 
should be noted that this method has a drawback - the information stored in registers are losing 
credibility. 


registers IFC 
Procedure of 


ESM sensor 
Sensor 2 and radar 

report fusion 

for three bba 


Sensor 1 


vectors 


Figure 3. The diagram of the information fusion process in the information fusion center IFC for two 
sensors and electronic entity record from IFC database. Explanations: mi —- BBA measure vector of i- 
th sensor, mr — generalized BBA measure vector that is a part of the electronic entity record in IFC, 
EER - electronic entity record in IFC database. 


In this case, the information fusion procedure performed in the IFC is carried out in accordance 
with the following formula: 


my = Rr(mp, m,,m,) y (7) 


wherein m; is a vector of the generalized BBA measure determined by the Ry rule based on 
the previous generalized BBA measure vector m,; and the new BBA measure vectors m, and m, 
sent by both sensors. The diagram of the identification information fusion from the ESM sensors is 
shown in Figure 3. 

In the further part of the paper, the combination rules of the BBA vector from the i - th sensor 
and the generalized BBA vector in the CFI are described. 


3.2. The rules of combination of BBA measures vectors 


This section presents formulas defining various combination rules for calculating basic belief 
assignments for the system shown in Figure 2 and Figure 3. The general forms are described in details 
in [4,11,12]. The information fusion rules of the DSm theory are presented below with the following 
constraints: 

e the properties of a set of hypotheses are described by the formulas (1) and (2), 

e for the first scheme (Figure 2), the information fusion procedure handles two information inputs: 
on one input, reports from two ESM sensors appear alternately, on the second input, electronic 
entity records from IFC database appear, 

e for the second scheme (Figure 3), the information fusion procedure handles three information 
inputs: on the first input, reports from a combined primary and secondary surveillance radar 
appear, on the second input, reports from a ESM sensor appear, on the third input, electronic 
entity records from IFC database appear. 


Dempster’s rule 


Dempster’s rule [13,14] of the BBA measure vector m; sent by the i-th sensor and the 
generalized BBA measure vector m, sin IFC is described for each 6; € © by the following formula: 
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l=... 
i = A 6xNO |= 0; 
ae a PR CST TC 
aed 
mp (8;)mi(9;) _ mpi(9) (8) 


~— 1-Y8_, Yo. me (Om (0)  1—key’ 
l#k 


wherein the k;; degree of conflict is defined by the formula: 


6 6 
kp = y. Mp (Ox)m; (81) = y y Mp (O,)m; (4) , 
k407.6 P= (9) 
1=1,...,6 l#k 
0 jNOK=0 
while 
One could notice that 
». mp (Ox)m;(0,) = 1, (11) 
k=1,...,6 
1=1,...,6 


k=T,...6 k=1,...6 
1=1,..6 1=1,..46 
6 jNOx=0 0jNOn#0 
6 6 6 6 
= YY: me(@i)mi(8.) + YY me(Eq)M4(G1) = 
k=11=1 k=1 l=1 
l#k l=k 
6 6 6 
= 6,.)m;(,) + 6,.)m; (Ox) = 1 (12) 
mp (8, )m; (81) Mp (8,)m (Ax) : 
k=1 [=1 k=1 
l#k 
From (12) it follows that if 
6 
Y, me Oe) U(e) = Lie. Y) mp()m4(6) = 0, (13) 
k=1 k=T,...6 
1=1,..,6 
then the degree of conflict is full. 
If 
6 
> Mp (8, )m; (Ox) = 0 , ie. y mp (8,)m,(6,) = 1, (14) 
k=1 k=1,...,6 
1=1,..,6 


then there is no conflict. 
mp(.) is Dempster-Shafer fusion result if and only if the denominator of the expression (8) is 


non-zero, i.e. the degree of conflict kp; is less than 1. 


The Proportional Conflict Redistribution rule PCR1 

PCR1 rule is the simplest and the easiest version of proportional conflict redistribution rule. The 
concept of the PCR1 rule assumes the calculation of the total conflicting mass (not worrying about 
the partial conflicting masses). The total conflicting mass is redistributed to all non-empty sets of 
hypotheses proportionally with respect to their corresponding non-empty column sum of the 
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associated mass matrix. The PCR1 rule is defined for every non-empty hypothesis in the following 


way: 
Cri cn(@) i 
mpr(6;) = Mpcri(9;) = ». mp(0x)m;(0,)] + — kp = 
k=1,..,6 
1=1,....6 
8x8 |=6; 
snip ene (15) 


= Mp (Gj); (O;) + — keri = Mp (j) + | Keri 


where c,;(6;) is the non-zero sum of the column saaneniig to the hypotheses 6; in the mass 


matrix 
M = laa (16) 
specified by the formula 
Cri(O;) = mMp(9;) + m;,(G;) (17) 


where: 

¢ mm, (i=1,2) is a row vector of the basic belief assignments masses of the i-th sensor's hypotheses , 
¢ mz, is arow vector of the basic belief assignments masses of the IFC system’s hypotheses, 

e =k; is the degree of mass conflict specified by the formula 


oe Mr (Ox)m;(O,) = . mr (Ox)m (81) , (18) 


peees k=11=1 
= ed : l#k 


e = dp;_ is the sum of all non-zero column sums of all non-empty sets 


6 


6 
i= > [me(6)) + mi(6)] = >. cri(B) (19) 


j=l j=1 


In our case dz; = 2 because 


6 6 
>, m2) =) mG) =1 (20) 
j=l j=l 
In addition 
Mri (6;) = Mp (8; )mM; (6;) (21) 


The Proportional Conflict Redistribution rule PCR2 


In PCR2 rule, the total conflicting mass k,; is distributed only to the non-empty sets involved 
in the conflict (not to all non-empty sets) and taken proportionally with respect to their corresponding 
non-empty column sum. 

A non-empty set 6; € @ is considered involved in the conflict if there exists another set 0; € 0 
which is neither included in @, nor includes 0, such that 6,9 6,=@ and m,;(0,.6,) > 0. The 


PCR2 rule is defined for every non-empty hypothesis 6; € @ in the following way: 


ma(@) = MoceatB) = Dy me(x)m(B)] + CG) EP ky = 
= cs i 
OKN81=6; 
= mp(8;)m;(8;) + C(8)) — - kp = Mpi(8;) + C(Gj)—— — kei (22) 


where 
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1, if 6; is involved in the conflict, 


otherwise. 


(23) 


Formula (23) can be written differently in the form (25), taking into account the definition of 


involvement in a conflict and formula (24) [12]: 


Mp (Bj NOx) = mp (8;) -m; (Ox) + Mp (Ox) - m;(8;) 


otherwise. 


1, if 4 Ox E 0,k #i] : mr (9;) ° m; (Ox) + Mp (Ox) * m;(9;) >0 


(24) 


(25) 


is the non-zero sum of the column corresponding to the hypotheses in the mass matrix M (16) 


specified by the formula 
Cri (Oj) = Mp(G;) + m,(G;) 


where: 


(26) 


¢ mm, (i=1,2) is a row vector of the basic belief assignments masses of the i-th sensor's hypotheses , 
¢ mz, is arow vector of the basic belief assignments masses of the IFC system’s hypotheses, 


e kg; is the degree of mass conflict specified by the formula (18), 


¢ eg; is the sum of all non-zero column sums of all non-empty sets only involved in the conflict 


6 6 
ei= >) [mr(G) +m(G)] =) cni(6)) = 2 C(8;)[m,(6,) + m,(6;)] = 2 C(6;)- (27) 


j€CF jECF 
where 

CF = {j =1,...,6: V0, € 0 mpi(6;NO,) > 0} 
and m,;(6; N 6) is defined by (24). 
In addition 


Mp (0;) = Mp (8;)m;(4;) 


(28) 


(29) 


It will be shown below that in the case of data used in numerical experiments (section 6) ep; = 2 
, this means that the PCR2 rule is equivalent to the PCR1 rule. The BBA vectors used there contain 


values less than 1, which means that 


(30) 


It follows that each BBA vector contains at least two non-zero components, that is Aj = 


1,...,6, 4k =1,...,6 with k #7 such that 


0 <mp(6;)<1 A 0< mp(O,)<1 


0<m,(8)<1 A 0<m(O,)<1 


(31) 


(32) 


From (31) and (32) it follows that if m,(6;)>0 , then there exists at least one value k # j such 


that m;(6,)>0 , which can be written in the following form 
Vj =1,...6: 0<mp(@)<1 > 3k=1,...,.6k#f: 0<m(%)<1 

From (33) it follows that 

Vj=1,...,.6: 0<m-(6)<1 => 3k =1,...,6,k #j: mp(Oj)-m;(Ox)>0 
The same applies: 

Vj=1,....6: 0<m(@)<1 > Ak=1,...,6,k Aj: m;(O;))+mp(O,)>0 
Taking into account (34), (35) and (25) one can obtain 

Vj =1,..,6: 0<m,(6;)<1 => 3k =1,...,6,k #j such that 


(33) 


(34) 


(35) 


(36) 
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Mp (O;) +m; (Ox) + m(G;) - Mp(Ox)>0 


Vj =1,..,6: 0<m(6;,)<1 > Ak =1,...,6,k #j such that 


m;(0;) + mp(Ox) + Mp (8;) +m; (Ox) >0 (37) 
From (36) and (37) it follows that 
Vj =1,...,6: 0<mp(6)<1 > C(6)=1 (38) 


This means that any hypothesis with a non-zero BBA value for any of the two sensors is involved 
in a conflict. 

From (27) it follows that 

6 6 6 
ei = > C(B)lmrr(B)) +mi(8))] =) C(B,)me(G) + > CG)mu(G}) (40) 
j=1 j=l j=l 

Using (36), (37) and (40), the value e;; will be determined. 

Because 


6 6 6 6 
DCO) G)= LY CO)meG) + DY CO)meG)= Yo me(G)=1 GD 
siete me(@,)=0 mn )>0 
6 6 6 6 

Cm (G)= YC) + YD C)m(G)= Y mG)=1 (42) 
=1 j=1 j=1 j=l 


J 


mj,(0;)>0 mj(@;)=0 mj;(0;)>0 
we get 
6 6 
emi =) C(B)me()) + > C(;)m(G}) = 2 (43) 
j=l j=l 


Considering (43), it can be said that in this case the PCR2 rule is equivalent to the PCR1 rule. For 
this reason, the results of the PCR2 rule are not presented in section 6, as they would be identical to 
the results of the PCR1 rule because we work only with Bayesian BBAs in this application. 


The Proportional Conflict Redistribution rule PCR3 


In PCR3 rule, one distributes the partial conflicting masses, instead of the total conflicting mass 
ky; , to the non-empty sets involved in the partial conflict. If an intersection is empty, for instance 
0,16, =@ , then the mass m(@;M @,) of the partial conflict is transferred to the non-empty sets 6; 
and @, proportionally with respect to the non-zero sum of masses assigned to 6, and respectively 
to 6, by the BBAs m,(.)and m;,(.). The PCR3 rule works if at least one set between 6, and @, is 
non-empty and its column sum is non-zero. 

The PCR3 rule is defined for every non-empty hypothesis 0; € @ in the following way: 


m()) = Moces() =E >. mr(x)miCOd] + [cre YSPC; O4)] = 
k=1,...,6 k 


=1,..., =1,...,6 
1=1,...,6 0,N8 ;=0 
O4N6{=8 ; 
= mpi(B) + [eni(8)) Y. SPE), 64)] (44) 
k=1,...,6 
k#j 


where 
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Mp (Ox)m;(6;) + Mp(G; mM; (Ax) 
SEER? (9;,0,%) = Cri (Oj) + Cri (Ax) 
0 for Cri (9;) + Cri(Ox) = 0 (45) 


for Cri(G;) + Cpi(O,) #0 


Cri (9; ) is the non-zero sum of the column corresponding to the hypotheses 6; in the mass matrix 
M (16) specified by the formula 


Cri(;) = mp (G;) + m,(G;) (46) 


The Proportional Conflict Redistribution rule PCR4 


The PCR4 rule redistributes the partial conflicting masses only to the sets involved in the partial 
conflict in proportion to the non-zero mass sum assigned to 6, and 6, by the conjunction rule 
according to the following formula: 


m,(6;) = Mecra(9;) = mp; (9;) + mp; (9;) bs Sr (6;, Ox) = 


= mpi(8j) + mei(8)_ > SEP**(6), 6x) (47) 


k=1,...,6 


where 
Mp (8; N Ox) 
Mp; (8;) + Mp; (Ox) 


for Cri (Gj) + Cp;(O,) # Oand Mp (9) Mp; (Ox) #0 


SPCR4(g. 6,) = i(0; 0 48 
FL (8, 6%) aD for Cpi(0;) + Cri(Ox) # 0 and mp;(6;) - Mp;(Ox) = 0 Os) 
0 for Cri(9;) + Cri(Ox) =0 
wherein 
Mri (9; NOx) = Mp (Ox)M; (Gj) + Mp (G;)Mi (Ox) (49) 
Mri (O;) = Mp (Bj) -m;(8;) (50) 
Cri (Gj) = Mp (Gj) + m,(G;) (51) 


If at least one of BBAs M,(. )or m;(. )is zero, the fraction is discarded and the mass Mp;(6; N 
6) is transferred to 6; and 8% proportionally with respect to their non-zero column sum of 


masses Cp; (6; ). 


The Proportional Conflict Redistribution rule PCR5 for two BBAs (two sources) 


Similarly to PCR2-PCR4 rules, PCR5 redistributes the partial conflicting mass to the hypothesis 
involved in the partial conflict. PCR5 provides the most mathematically precise [4,11,12] 
redistribution of conflicting mass to non-empty sets in accordance with the logic of the conjunctive 
rule. However, it is more difficult to implementation. The PCR5 rule is defined for every non-empty 
hypothesis 0; € @ in the following way: 


mic(8,) = Mpces (Bj) = Mer()) + > SEF (6,64) = mri(Gj) + > SEPR(G;,6,) (52) 
= = 6 


where 
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mp (6;)" -mj(8x) m;(9;)? > Mp (Ox) ive 

m,(6;) +m;(6,) m;(6;) + mp (Ox) (53) 
Mp (6;) +m;(0,) #0 and m;(6;) +mr(0,) #0 

0 for mp (8;) +m;,(0,) = 0or m;(8;) +m,p(O,) = 0 


Sit”? (6), 0%) = 


wherein 
Mri (0;) = mp(8;) -m;(9;) (54) 


In the formula (52), the component S7£*° is equal to zero if both denominators are equal to zero. 
In the formula (53), if a denominator is gate then component is discarded. 


The Proportional Conflict Redistribution rules PCR5 and PCR6 for three BBAs (three sources) 

In [11,12] improved proportional conflict redistribution rules of combination of basic belief 
assignments PCR6, PCR5+ and PCR6+ are presented. The authors point out that these rules should 
be applied if and only if we are to combine more than two BBAs. If we have only two BBAs to combine 
(s = 2) we always get mecrs = mpcrs+ = mpcro = mecro+ because in this case the PCR5, PCR5+, PCR6, and 
PCR6+ rules coincide. Below are the formulas that define the PCR5 and PCR6 rules for 3 BBAs. 

The PCR5 rule for three BBAs (three sources) is defined for every non-empty hypothesis in the 


following way: 


mp(8;) = Moca (8) = a (55) 
FL) = Mpcrs\4y) = Fe m'(6;) 
t=1 
wherein 
(Bj) = mer2(8j) + >. SEER} 0,81) + > SLPESS(6;, 6x) + 
k=1,...,6 k=1,...,6 
1=1,...,6 04.8 ;=0 
818 {NG ; =0 
+ S2F12° (Bj, Ox) = 
k=1,...,6 
O46 ;=0 
= me12(8)) + YY SEF (6,844) + ps POSS (6;,8,.) + 
k=1,..,6 1=1,..,6 k=1.,..., 
k#j l#jAltk po 
+ )° S2EG8 (6,4) = 
k=1,...,6 
jtk 
= me12()) +). | >. SEEFS(6), 64,61) + STESES(6;, Ox) + S2ES8S(6;,O%)] (66) 


k#j Liejale#k 


me (Bj)? «m4 (H.) + M21), M(H.) = M(H)? + M2 Gx) 
me(8;) + m, (6) + m2(8.) me(H.) +m, 6) + mG.) 

me (By.) My (6) + M2 (8)? (57) 
mie (By) +m, (6,) + m3 (6) 


She? (6;, O,,8;) = 


pens (a, g,) = ACN MaCOu)-ma(Gx) , me(Ox) “ms (G;)? ma(Or) 

trai Ped) = TG) +m (Gy) + Mz Oq) | MeO) + 4G) + m2 Oy) 
me(j.) “4 (84.) “13(6;)? (58) 
mie (By.) + my (By) + m8) 
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S2PSRS (6, 6,) = mp (8;)* +m, (6;)? -mM2(Ox) | M(x) -™71(8))? + m2(;)? 
“< mp(9;) + m,(9;) + m2 (A) Mp (A) +m, (9;) + m2 (4;) 
mp (8j)? +m (Ox) + m2 (8;)* (59) 
mp (9;) +m, (Ox) + m2(4;) 
ee _ mG) 
Mp (9;) = Mpcrs (Gj) = yo, m’(6,) (60) 


In the formulas (57)-(59), if a denominator is zero, then component is discarded. 
The quotient in the formula (55) ensures the normalization of the BBA vector m;, which ensures 


6 6 
>»: mp (6;) = y Mpcrs(9;) = 1 
i=1 i=1 


The PCR6 rule for three BBAs (three sources) is defined for every non-empty hypothesis 6; € 0 
in the following way: 


that 


m (6;) 


my(0;) = Mpcps (8) = =~ (61) 
FY] PCRS\Yj em (0;) 
wherein 
m'(6;) = Mpr2(6;) + >). SEER(G}, 0,61) +). STBESS(6;,6,) + 
k=1,..., : k=1,...,6 
[=i 640 j=0 
xn 08, =@ 
a S2643° 8}, 4) = 
k=1,...,6 
O68 ;=9 
= mer2(8)+ YD. SEFC). 8x, 81) + 2 S158° (8,4) + 
k=1,..,6 1=1,...,6 k=1,.., 
k#j ltjaltk ng 
+ D) S27SR(6,,6,) = 
k=1,...,6 
j#k 
= Mp12(9;) + y » SEIS? (Bj, Ox 1) + S1PSF°(8;, Ox) + S2EEF° (Gj, Oi.) (62) 
k=1,..,6] 1=1,...,6 
k#j Liejaltk 
with 


mp(8j)? +m (Bie) M2 (81) Me(1) +4 (8))? Mz (Gx) 


PCRS = m.(0,) +m.(0.) + m.(6,) 
SPIE (9, Pi 1) = mp(8;) +m, (Ox) + m2(4,) mp(4,) +m, (4;) + mz (Ox) 


mye (Ox), (8,) -m2(8,)? (63) 
mie (Ox) + m4 (81) + (8) 


mp(8j)? + m1 (Bie) M2) M(x) + M1 (Bj) + M2 (Ax) 


5 19085(6,,0,) = MGV” Mao) ma Oo , mnOr) mG mal) 
Mp (6;) +m, (Ox) +™2(O%) mp (O,) + m,(6;) +m (6x) 


mg (Ox) M4 (8y.) + M9 (8;)? (64) 
mie (Ox) + 3 (Ox) + m2 (4)) 


mp (8)? «m4 (8) + M2 (Ox) + Me(4j) + mM (6;)* M2 (x) a (65) 
ei) May Tea) 


crn ran? x rns x cnn? 


S2F 12° (8), 9%) = 
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Peli (8;)? my (8;.) «M2 (B;) + M(B) +My (Ox) + M2 (8;)” 
mp (Oj) + mM, (Ox) + m2 (G;) 


and 
™Mp12(9;) = Mp(G;) -m,(8;) - m2(4;) (66) 


In the formulas (63)-(65), if a denominator is zero, then component is discarded. 
The quotient in the formula (61) ensures the normalization of the BBA vector m; , which 


ensures that 
6 6 
>. me() = >. Meces(6,) = 1 
i=1 i=1 


Comparing the two fusion schemes (Figures 2 and 3), it should be noted that sequential and 
global information fusion generally produces different results [4], i.e. 


PCRS(mp,™m,,M2) # PCRS5(PCR5(mp, m,), M2) # PCRS(PCR5(mp, m2), m,) (67) 


In addition, the article experimentally verified the theorem on the inequality of the results of 
both PCR5 and PCR6 rules for three BBAs (three sources) presented in [4]: 


PCRS(mgz, m1,Mz2) # PCR6(Mp, mM,,M2) (68) 


4. Basic belief assignment for combined primary and secondary surveillance radars 


The paper assumes that the analyzed radar sensor consists of two radars: primary and 
secondary. Therefore, the probability of correct detection and correct identification of a target is 
expressed by the following formula: 


Poi = Pa’ Pirg (69) 


where Pd is the probability of correct detection of target by a primary radar and PIFF is the probability 
of correct reply for interrogation. If a target has been detected by the primary radar and there is a lack 
of proper identification by the secondary radar one can assume that the target has a value of attribute 
identification of UNKNOWN - U. So, one can write the following relation: 


m(U) = Pg:(1—- Pree), (70) 


where m(U) is the mass of probability for a value of UNKNOWN identification attribute. 

A method of calculating the probabilities Py and Piz, is presented in [5,15,16]. 

The way of allocation of the remaining mass of probability (1 - m(U)) will be described in this 
section. One assumes there that every simulated target should have a base value of attribute 
identification from the set 


Ley ={N;.F,H,} (71) 


where: 

e Ns-—base NEUTRAL identity, 
e Fs—base FRIEND identity, 

e Hs—base HOSTILE identity. 

STANAG 1241 introduces in addition to the basic set of attribute identification values also 
secondary (additional) attribute identification values: SUSPECT (S) and ASSUMED FRIEND (AF). 
According to Figure 1 one can introduce a table of possible attribute values transitions between set 
(10) and the set of secondary attribute identification values: 


Zs = {Ns, Fs, Hs, AF, S} (72) 


Table 1. Transformation of the base belief assignment mass into the secondary belief assignment 
mass 
Base identification 
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Belief mass values contained in the Table 1 determine how the mass of the base belief assignment 
has been transformed into the mass of the secondary belief assignment. They can be estimated as 
empirical frequencies based on recorded archive events. 


Table 1. Transformation of the base belief assignment mass into the secondary belief assignment 


mass. 
Base identification > Fp Ne He 
Fs m(Fs| Fg) 0 0 
Ns 0 m(Ns|Np) 0 
Hs 0 0 m(Hs|Hg) 
AE m(AF|Fa) m(AF|Nz) 0 
S 0 m(S|Fp) m(S|Hp) 


Of course, have the normalization conditions satisfied: Yyez,,m(x|Fs)=1, Yxez,,m(x|Ng) = 1 
and )ixez,,m(x|Hg) = 1. 
The final values of the belief mass of secondary attribute identification values are calculated 
according to the formulas: 
1. For a target with the FRIEND base value of an attribute identification 
m(U) = Pal — Pree), 
m(AF) = m(AF|Fg)(1 — m(U)), 
m(Fs) = m(Fs|Fp)(1 -_ m(U)). 
2. For a target with the NEUTRAL base value of an attribute identification 
m(U) = Pal — Pree), 
m(AF) = m(AF|Ng)(1 — m(U)), 
m(S) = m(S|Ng)(1 —m(U)), 
m(Ns) = (1 — m(AFINg) — m(S|Ng))(1 — m(U)). 
3. For a target with the HOSTILE base value of an attribute identification 
m(U) = Pa(1 — Pree), 
m(Hs) = m(Hs|Hg)(1 al m(U)), 
m(S) = m(S|H,)(1 — m(U)). 
Other final values of the belief mass of secondary attribute identification values are equal zero. 


5. Basic belief assignment for ESM sensors 


An ESM sensor is a passive sensor that captures incoming electromagnetic signals generated 
first of all by radar emitters mounted on air or maritime platforms. This sensor recognizes radar 
signals determining values of their distinctive features. In this paper we will not deal with methods 
of radar signals recognizing in details. However, we will use information about these methods to 
identify platforms generating the signals according to STANAG 1241 — NATO Standardization 
Agreement and Dezert-Smarandache theory. As previously it was stated, we are interested in three 
basic values of identification: friend, hostile and neutral, and two secondary values: suspicious and 
assumed friendly. In addition, we will assume that in some situations it is not possible to determine 
the identity of the emitter carrier platform. To clarify this issue, we should briefly describe the method 
of determining the identification of the emitter carrier platform that generated the captured signal. 
The sensor recognition system is equipped with a database that can be divided into three 
components: a platform database, an emitter list and a geopolitical list [10]. The platform database 
(PDB) contains information about platforms that can be met in the area of interest along with their 
equipment with emitters. The emitter name list (ENL) includes all emitters corresponding to each 
platform of the PDB and contains the values of the signal distinctive features for each emitter. The 
values of distinctive features are the basis for the procedure of recognizing a captured signal. The 
geopolitical list (GPL) provides the allegiance of various countries and platforms and allows to 
identify them in accordance with STANAG 1241. 

The algorithm of signal recognition is realized in two stages: 
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1. Verification at the level of signal quality features. The second stage is executed after a positive 
assessment of the conformity of quality features. 

2. The signal recognition procedure determines the distances between the distinctive features of 
the recognized signal and the distinctive features of all pattern signals stored within the emitter list. 

Let us introduce the following notation: 

x, — vector of distinctive features of the recognized signal, 

x; — vector of distinctive features of i-th pattern signal (i— the number of the pattern signal, i € 
{1,...,M} ), 

d,; = d(xs,x;,) — the distance between the distinctive features vector of the recognized signal 
and the distinctive features vector of i-th pattern signal; the distance d,; is the Mahalanobis distance 
taking into account the correlations of the distinctive features. 

The signal recognition classifier compares the distance d(x,,x;) with the acceptable positive 
distance of the classification 6 . The distance 6 is the limit that we will interpret as a boundary of 
emitter pattern recognition. We will divide the set of pattern signals into two subsets: the patterns 
satisfying the positive classification condition in relation to the recognized signal s — D3 and the 
patterns that do not satisfy the positive classification condition - D;. The formal definition is as 
follows: 


Di ={ ie (1,...,M}ld,; < 6} (73) 


D; ={ie {(1,...,M}d,; > 5} (74) 


In the paper we propose the following method of determining the basic belief assignment on a 
set of pattern signals, which is related to the distance between a signal and a pattern in the distinctive 
features space: 


m@i=e (75) 
As one can see from the formula (75) if d(x,,x;) = 0 then m,(i) = 1 , whereas if d(x,,x;) > 0 
then 0 <_m,(i) <1 . The above measure is not normalized, hence we will normalize it 
ms (i) 
jai m,(i) 
The sum of the measures assigned to all the emitters, whose distinctive features lie outside the 
limit 6 , will be treated as a measure assigned to the base hypothesis “unknown” (U) 


Ms(t) = (76) 


fisU) = ) i. 7) 
ieD; 

To determine the belief measure of other base hypotheses (H, F, N) and secondary hypotheses 
(AF and S), we should introduce formal definitions of sets contained in the sensor database and used 
for recognition of captured signals. As it was mentioned above the set of all the necessary data for 
platform identification can be divided into three sets: PDB — a platform database, ENL — an emitter 
name list and GPL — a geopolitical list: 

PDB - the platform database contains information about all platforms observed in the area of 
interest, including information on all emitters mounted on each platform; we assume that one 
platform can have many emitters and the same type of emitters can be installed on many platforms; 
the PDB contains also information on the national affiliation of each platform, 

ENL -— the emitter name list is a set of information about all recognized emitters in the area of 
interest; this set contains the mean values of the distinctive features of emitter signals (so-called signal 
patterns) and their standard deviations, 

GPL - the geopolitical list contains base values of identification attributes (H, F, N) assigned to 
the various countries. 

We will also introduce additional notations used in this paper: 

PDBL - the list of platform numbers that are stored in the PDB, 

PL(i) — the set of numbers of platforms which have the emitter with number 


wo 
1 
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IPL(j) — the base identification attribute of the platform with number “j” determined on the basis 
of the information contained in PDB and ENL (IPL(j) € { F A, M ). 

The set of signal patterns satisfying the positive classification condition in relation to the 
recognized signal s denoted as DJ can be divided into disjunctive subsets according to the values of 
the carrier platform identification features: 


R= DP up ype opi’ ups, (78) 


Dit npt=0, k #1, k,l€ {FANAFS}. (79) 


Each subset of the set Dt for the base identification is defined as follows: 


DiF ={ie Dt|vj € PL() IPL() = F}, (80) 
Di# ={ie D?|vj € PL(i) IPL() = H}, (81) 
Di ={ie Dt |vj € PL(i) IPL(j) = N}. (82) 


In a similar way, one can define subsets of the set D} for the secondary identification (AF, S): 


D4? ={ie Dt|aj € PL(i) IPL(j) =F A 4j € PL(i) IPL) =N} (83) 


Dis ={ie Dt:4j € PL(i) IPL(j) =H A 4j € PL(i) IPL) =N} (84) 


One can notice we assume in this paper that no emitter type can be installed simultaneously on 
platforms with identifications F and H: 


{ie Dt|aj¢€PL(i) IPL(j)=F A Aj€PL(i) IPL(j)=H}=0 (85) 


Introducing the definition of subsets of the set allows us to determine the belief masses for all 
identification features: 


m,(F) = oH m,(i) , m,(H) a 2 m,(i) , m,(N) = 2 m,(i) , (86) 
m,(AF) = Lm m,(S)= 2 m, (i) : (87) 


It should be emphasized that the method presented here is different than in [6,17]. These papers 
assume that ESM sensors can only generate basic declarations with attribute values FRIEND, 
HOSTILE and NEUTRAL but in this paper, we assume, that ESM sensors can generate declarations 
from an extended set of attribute values (additionally ASSUME FRIEND, SUSPECT and 
UNKNOWN). 


6. Numerical experiments of fusion of identification information from ESM sensors 


6.1. Simulation scenarios 


The paper [6] presents a typical simulation scenario for testing the identification information 
fusion. The authors formulated several requirements that should be met by sucha scenario. It should: 
1) adequately represent the known ground truth of the emitter identification, 

2) include sufficient numbers of incorrect associations to be realistic and to test the robustness of 
the rules to temporary incorrect sensor decisions, 

3) provide only partial knowledge about the ESM sensor declarations, and thus contain 
uncertainty, 

4) allow to show stability in case of countermeasures, 

5) allow to switch identification when the ground truth changes. 

The authors [6] propose the following parameters of the scenario: 

1) ground truth of identification is FRIEND (F) for the first 50 iterations of the scenario and 

HOSTILE (A) for the last 50 iterations, 
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2) the number of correct associations is 80% of all iterations, the number of incorrect associations 
caused by countermeasures is 20% of all iterations in a randomly selected moments of time, 
3) ESM sensor declarations have a mass of 0.7 for the most credible identification and 0.3 for the 

identification of UNKNOWN (U). 

The assumption 5) is not considered in this paper, assuming that the real object does not change 
its real identity while performing the mission. Therefore, assumption 1) regarding scenario 
parameters becomes obsolete. 

The following assumptions concerning the parameters of the scenario are made in this paper: 
1) the real value of identification is constant in each scenario and is equal to FRIEND (F) — in the 

scenarios 1, 2 and 5 and HOSTILE (A) - in the scenarios 3, 4 and 6; 

2) the above declarations are transmitted by sensor number 1 with the real identification mass 

equal to 0,7 and the mass of complementary identification (UNKNOWN) equal to 0,3; 

3) the second sensor shall transmit its declarations in accordance with the Tables 2 and 3 for the 

scenarios 1 and 2 and with the Tables 4 and 5 for the scenarios 3 and 4. 


Table 2. Belief mass values for the second sensor for the scenarios 1 and 5. 


Type of identification F N H AF S Uu 
Correct identification (80% of events) 0.6 0.1 0 0.2 0 0.1 
Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 
Table 3. Belief mass values for the second sensor for the scenario 2. 
Type of identification F N H AF S U 
Correct identification (80% of events) 0.7 0.1 0 0.1 0 0.1 
Incorrect identification (20% of events) 0 0.1 0.7 0 0.7 0.1 
Table 4. Belief mass values for the second sensor for the scenario 3. 
Type of identification F N H AF S U 
Correct identification (80% of events) 0 0.1 0.6 0 0.2 0.1 
Incorrect identification (20% of events) 0.6 0.1 0 0.2 0 0.1 
Table 5. Belief mass values for the second sensor for the scenarios 4 and 6. 
Type of identification F N H AF S Uu 
Correct identification (80% of events) 0 0.1 0.7 0 0.1 0.1 
Incorrect identification (20% of events) 0.7 0.1 0 0.1 0 0.1 


One should note that scenario 2 differs from scenario 1 with a greater belief mass assigned to 
incorrect identification of the recognized emitter. The scenarios 3 and 4 are similarly different. 

Scenarios 1-6 for the sensor 1 have been presented in Figures 4 and 5. All deterministic scenarios 
for the sensor 2 have been presented in Figures 6-9. 

This section may be divided by subheadings. It should provide a concise and precise description 
of the experimental results, their interpretation, as well as the experimental conclusions that can be 
drawn. 
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Figure 4. The course of scenarios number 1, 2 and 5 for 


sensor 1. 
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Figure 6. The course of scenarios number 1 for sensor 2. 
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Figure 8. The course of scenarios number 3 for sensor 2. 
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Figure 7. The course of scenarios number 2 for sensor 2. 
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Figure 9. The course of scenarios number 4 for sensor 2. 


The Monte Carlo method of generating the scenario for the sensor 2 is also used in this paper. 
Moments in which incorrect identifications occurred are generate by the pseudorandom integer 
number generator from the range [0, 100]. Examples of scenarios are shown in Figures 10 and 11. 
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Figure 10. The course of Monte Carlo scenarios number 5 Figure 11. The course of Monte Carlo scenarios number 6 
for sensor 2. 
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6.2. Calculation results for deterministic scenarios 


Dempster’s rule 


Dempster’s rule is not resistant to a situation when the degree of conflict kri = 
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1. This means the 


total conflict between the mass vector sent by the sensor and the mass vector of the information fusion 
center, which occurs when each non-zero belief mass value sent by the sensor corresponds to zero 


belief mass value of the vector determined by the information fusion center and vice versa. 


The simulation results of the identification information fusion using Dempster’s rule have been 


presented for the deterministic scenarios 1 and 3 in Figures 12 and 13. 
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Figure 12. The values of the resulting belief mass for 
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scenario 1 and Dempster’s rule. 


The PCR1 rule 
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Figure 13. The values of the resulting belief mass for 
scenario 3 and Dempster’s rule. 


The simulation results of the identification information fusion using the PCR1 rule for the 
deterministic scenarios 1, 2, 3 and 4 are presented in Figures 14-17. 


mass value m(.) 


0 20 40 60 80 100 120 140 160 180 200 
time step 


Figure 14. The values of the resulting belief mass for 
scenario 1 and the PCRI rule. 
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Figure 16. The values of the resulting belief mass for 
scenario 3 and the PCRI rule. 
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Figure 15. The values of the resulting belief mass for 
scenario 2 and the PCRI rule. 
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Figure 17. The values of the resulting belief mass for 
scenario 4 and the PCRI rule. 


The simulation results of the identification information fusion using the PCR3 rule for the 
deterministic scenarios 1, 2, 3 and 4 are presented in Figures 18-21. 
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Figure 18. The values of the resulting belief mass for 
scenario 1 and the PCR3 rule. 
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Figure 19. The values of the resulting belief mass for 
scenario 2 and the PCR3 rule. 
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Figure 20. The values of the resulting belief mass for 
scenario 3 and the PCR3 rule. 
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Figure 21. The values of the resulting belief mass for 
scenario 4 and the PCR3 rule. 
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The simulation results of the identification information fusion using the PCR4 rule for the 
deterministic scenarios 1, 2, 3 and 4 are presented in Figures 22-25. 
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Figure 22. The values of the resulting belief mass for 
scenario 1 and the PCR4 rule. 
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Figure 24. The values of the resulting belief mass for 


scenario 3 and the PCR4 rule. 
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Figure 23. The values of the resulting belief mass for 
scenario 2 and the PCR4 rule. 
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Figure 25. The values of the resulting belief mass for 
scenario 4 and the PCR4 rule. 


The simulation results of the identification information fusion using the PCR5 rule for 2 BBAs 
for the deterministic scenarios 1, 2, 3 and 4 are presented in Figures 26-29. 
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Figure 26. The values of the resulting belief mass for 


scenario 1 and the PCR5 rule for 2 BBAs. 
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Figure 28. The values of the resulting belief mass for 
scenario 3 and the PCR5 rule for 2 BBAs. 
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Figure 27. The values of the resulting belief mass for 
scenario 2 and the PCR5 rule for 2 BBAs. 
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Figure 29. The values of the resulting belief mass for 
scenario 4 and the PCR5 rule for 2 BBAs. 


The simulation results of the identification information fusion using the PCR5 rule for 3 BBAs 
for the deterministic scenarios 1, 2,3 and 4 are presented in Figures 30-33. 
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Figure 30. The values of the resulting belief mass for 
scenario 1 and the PCR5 rule for 3 BBAs. 
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Figure 31. The values of the resulting belief mass for 
scenario 2 and the PCR5 rule for 3 BBAs. 
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Figure 32. The values of the resulting belief mass for 


scenario 3 and the PCR5 rule for 3 BBAs. 
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Figure 33. The values of the resulting belief mass for 
scenario 4 and the PCR5 rule for 3 BBAs. 


The simulation results of the identification information fusion using the PCR6 rule for 3 BBAs 
for the deterministic scenarios 1, 2,3 and 4 are presented in Figures 34-37. 


mass value m(.) 


TTT TT 


ee mr a eR A | | 


0 20 40 60 80 100 120 140 160 180 200 
time step 


Figure 34. The values of the resulting belief mass for 
scenario 1 and the PCR6 rule for 3 BBAs. 
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Figure 36. The values of the resulting belief mass for 
scenario 3 and the PCR6 rule for 3 BBAs. 
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Figure 35. The values of the resulting belief mass for 
scenario 2 and the PCR6 rule for 3 BBAs. 
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Figure 37. The values of the resulting belief mass for 
scenario 4 and the PCR6 rule for 3 BBAs. 


The presented results (Figures 12-37) allow us to conclude that the applied methods of managing 
conflicts in the information fusion allow to draw correct conclusions about the real identification of 


the recognized object. 
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The application of the decision threshold for the belief mass at the level mo=0,37 for the PCR1 
rule and ma =0.45 for the PCR3, PCR4 and PCR5 rules for scenarios 1 and 3 allows us to properly 
evaluate the identification of the recognized object: for scenario 1 — FRIEND and for scenario 3 — 
HOSTILE. For scenarios 2 and 4, the optimal thresholds are ma =0.4 for the PCR1 rule and ma =0.45 
for the PCR3, PCR4 and PCR5 rules respectively. When assessing the interval between the minimum 
resultant mass for correct identification and the maximum resultant mass for misidentification, the 
worst results are reached by the PCR1 rule and the rules of PCR3, PCR4 and PCR5 behave similarly 
and are better than the rule PCR1. 

The research carried out for deterministic scenarios shows that the PCR5 rule for 3 BBAs and the 
PCR6 rule for 3 BBAs behave very similarly. They restores the correct identification after the 
occurrence of temporary misidentification much faster than the rules PCR1 — PCR5 for 2 BBAs 


6.3. Calculation results for Monte Carlo scenarios 


Dempster’s rule 


In the Monte Carlo scenario, Dempster’s rule behaves similarly to a deterministic scenario. It is 
not resistant to a situation when the degree of conflict kFi = 1. This means the total conflict between 
the mass vector sent by the sensor and the mass vector of the information fusion center which occurs 
when each non-zero belief mass value sent by the sensor corresponds to zero belief mass value of the 
vector determined by the information fusion center and vice versa. 

The simulation results of the identification information fusion using Dempster’s rule are 
presented for the Monte Carlo scenarios 5 and 6 in Figures 38 and 39. 

rule have been presented for the deterministic scenarios 1 and 3 in Figures 12 and 13. 
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Figure 38. The values of the resulting belief mass for Figure 39. The values of the resulting belief mass for 
Monte Carlo scenario 5 and Dempster’s rule. Monte Carlo scenario 6 and Dempster’s rule. 
The PCR1 rule 


The simulation results of the identification information fusion using the PCR1 rule for the Monte 
Carlo scenarios 5 and 6 are presented in Figures 40 and 41. 
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Figure 40. The values of the resulting belief mass for Figure 41. The values of the resulting belief mass for 
Monte Carlo scenario 5 and the PCR! rule. Monte Carlo scenario 6 and the PCR1 rule. 
The PCR3 rule 


The simulation results of the identification information fusion using the PCR3 rule for the Monte 
Carlo scenarios 5 and 6 are presented in Figures 42 and 43. 
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Figure 42. The values of the resulting belief mass for Figure 43. The values of the resulting belief mass for 
Monte Carlo scenario 5 and the PCR3 rule. Monte Carlo scenario 6 and the PCR3 rule. 
The PCR4 rule 


The simulation results of the identification information fusion using the PCR4 rule for the Monte 
Carlo scenarios 5 and 6 are presented in Figures 44 and 45. 
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Figure 44. The values of the resulting belief mass for Figure 45. The values of the resulting belief mass for 
Monte Carlo scenario 5 and the PCR4 rule. Monte Carlo scenario 6 and the PCR4 rule. 
The PCR5 rule for 2 BBAs 


The simulation results of the identification information fusion using the PCR5 rule for 2 BBAs 
for the Monte Carlo scenarios 5 and 6 are presented in Figures 46 and 47. 
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Figure 46. The values of the resulting belief mass for Figure 47. The values of the resulting belief mass for 


Monte Carlo scenario 5 and the PCR5 rule for 2 Monte Carlo scenario 6 and the PCR5 rule for 2 
BBAs. BBAs. 
The PCRS5 rule for 3 BBAs 


The simulation results of the identification information fusion using the PCR5 rule for 3 BBAs 
for the Monte Carlo scenarios 5 and 6 are presented in Figures 48 and 49. 
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Figure 48. The values of the resulting belief mass for Figure 49. The values of the resulting belief mass for 


Monte Carlo scenario 5 and the PCR5 rule for 3 Monte Carlo scenario 6 and the PCR5 rule for 3 
BBAs. BBAs. 
The PCR6 rule for 3 BBAs 


The simulation results of the identification information fusion using the PCR6 rule for 3 BBAs 
for the Monte Carlo scenarios 5 and 6 are presented in Figures 50-51. 
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Figure 50. The values of the resulting belief mass for Figure 51. The values of the resulting belief mass for 


Monte Carlo scenario 5 and the PCR6 rule for 3 Monte Carlo scenario 6 and the PCR6 rule for 3 
BBAs. BBAs. 


The presented results show that due to the high intensity of sending reports with incorrect 
identifications in the middle part of the scenarios, the information fusion rules (apart from the 
PCR4,PCR 5 and PCR6 rules) determine the maximum resulting mass for incorrect identification. The 
PCR5 for 3 BBAs and PCR6 for 3 BBAs rules are the fastest to restore the correct identification after 
receiving several incorrect reports. 


7. Numerical experiments of fusion of identification information from ESM sensors and radars 


7.1. Numerical experiments scenarios 


We assume that we will combine attribute information from 2 sensors: a combined primary and 
secondary surveillance radar and an ESM sensor. These sensors work asynchronously. Upon receipt 
of the sensor's declaration in the form of a vector of masses, we fuse this vector with the vector of the 
actual values of the declaration masses for the fuser’s frame of discernment. The frequency of 
transmission of sensor declarations depends on the rules of the data exchange network and on the 


doi:10.20944/preprints202306.1892.v1 


28 


technical characteristics of the sensors. Various combination methods are presented in [13, 14]. In this 
paper, two of the methods of proportional redistribution conflict (PRC5 and PCR6 [11,12]) has been 
used. The numerical model of combined primary and secondary surveillance radars was taken from 
[15,16]. 

Numerical experiments have been performed for the following data: 

for combined primary and secondary surveillance radars sensor: 

Pra=10~°, Rmax= 100 [km], P7=0.7, o¢=2 [m2], Prpp= 0.962 

and the following table of masses (compare Table 1): 


Table 6. Transformation of the base belief assignment mass into the secondary belief assignment mass 
for combined primary and secondary surveillance radar. 


(Scenario Nr, Base identification) > (1,Fs) (2,NB) (3,Hs) 
Fs 0.8 0 0 
Ns 0 0.5 0 
Hs 0 0 0.7 
AF 0.2 0.3 0 
S 0 0.2 0.3 


The flight path of air object was 30 km away from the sensor (in the horizontal plane), flight 
altitude was 1 km and radar cross-section was 1 sq. m. 
e for EMS sensor: 
the following assumptions concerning the parameters of the scenario for ESM sensor are made 
in this paper: 
1) the real value of identification is constant in each scenario and is equal to FRIEND (F) — in the first 
scenario and HOSTILE (H) - in the second scenario; 
2) the above declarations are transmitted by sensor number 1 with the real identification mass equal 
to 0,7 and the mass of complementary identification (UNKNOWN) equal to 0,3; 
3) the second sensor shall transmit its declarations in accordance with the Tables 1 and 2 for the 
scenarios 1 and 2 and with the Tables 3 and 4 for the scenarios 3 and 4. 
In the Tables 7-12 the mass values for all possible declarations for six scenarios for ESM sensor 
are presented. 


Table 7. Belief mass values for the second sensor (ESM) for the scenario 1. 


Type of identification F N H AF S Uu 
Correct identification (80% of events) 0.6 0.1 0 0.2 0 0.1 
Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 


Table 8. Belief mass values for the second sensor (ESM) for the scenario 2. 


Type of identification F N H AF S Uu 
Correct identification (80% of events) 0 0.5 0.3 0 0.2 0 
Incorrect identification (20% of events) 0 0.4 0.2 0 0.3 0.1 


Table 9. Belief mass values for the second sensor (ESM) for the scenario 3. 


Type of identification F N H AF S Uu 
Correct identification (80% of events) 0 0.1 0.7 0 0.1 0.1 
Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 


Table 10. Belief mass values for the second sensor (ESM) for the scenario 4. 


Type of identification F N H AF S 
Correct identification (80% of events) 0.1 0.7 0.1 0 0 0.1 
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Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 


Table 11. Belief mass values for the second sensor (ESM) for the scenario 5. 


Type of identification F N H AF S Uu 
Correct identification (80% of events) 0.6 0.1 0 0.2 0 0.1 
Incorrect identification (20% of events) 0 0.1 0.6 0 0.2 0.1 


Table 12. Belief mass values for the second sensor (ESM) for the scenario 6. 


Type of identification F N H AF S u 
Correct identification (80% of events) 0.1 0.7 0.1 0 0 0.1 
Incorrect identification (20% of events) 0.6 0.1 0 0.2 0 0.1 


Scenarios 1-6 for the sensor 1 have been presented in Figures 52,53 and 54. All deterministic 
scenarios for the sensor 2 have been presented in Figures 55-60. 
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Figure 52. The course of scenarios number land 4 for Figure 53. The course of scenarios number 2and 5 
sensor 1. for sensor 1. 
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Figure 54. The course of scenarios number 3and 6 for 
sensor 1. 
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Figure 55. The course of scenarios number 1 for 


sensor 2. 
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Figure 57. The course of scenarios number 3 for 


sensor 2. 
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Figure 59. The course of scenarios number 5 for 


sensor 2. 
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Figure 56. The course of scenarios number 2 for 
sensor 2. 
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Figure 58. The course of scenarios number 4 for 
sensor 2. 
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Figure 60. The course of scenarios number 6 for 
sensor 2. 


Scenarios 1-3 assume relatively small changes in the mass of all declarations. Scenarios 1-3 
assume significant changes in the credibility mass of all declarations (small errors). Scenarios 4-6 
assume significant changes in the mass of all declarations (large errors). 
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7.2. Calculation results for four proportional conflict redistribution rules 


The PCR5 rule for 2 BBAs 


The simulation results of the identification information fusion using the PCR5 rule for 2 BBAs 
for the deterministic scenarios 1-6 are presented in Figures 61-66. 
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Figure 61. The values of the resulting belief mass for 


Figure 62. The values of the resulting belief mass for 
scenario 1 and the PCRS rule for 2 BBAs. 
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Figure 63. The values of the resulting belief mass for 
scenario 3 and the PCR5 rule for 2 BBAs. 


Figure 64. The values of the resulting belief mass for 
scenario 4 and the PCR5 rule for 2 BBAs. 
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Figure 65. The values of the resulting belief mass for 
scenario 5 and the PCR5 rule for 2 BBAs. 
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Figure 66. The values of the resulting belief mass for 
scenario 6 and the PCR5 rule for 2 BBAs. 


For the PCR5 rule for 2 BBAs, the application of the decision thresholds at the belief mass level 
mo=0,43 for scenario 4, ma=0.35 for scenario 5 and muo=0. 5 for scenario 6 allows us to properly evaluate 
the identification of the recognized object for most time moments. 


The PCR5 rule for 3 BBAs 


The simulation results of the identification information fusion using the PCR5 rule for 3 BBAs 
for the deterministic scenarios 1-6 are presented in Figures 67-72. 
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Figure 67. The values of the resulting belief mass for 
scenario 1 and the PCRS rule for 3 BBAs. 
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Figure 69. The values of the resulting belief mass for 
scenario 3 and the PCR5 rule for 3 BBAs. 
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Figure 71. The values of the resulting belief mass for 
scenario 5 and the PCRS rule for 3 BBAs. 
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Figure 68. The values of the resulting belief mass for 
scenario 2 and the PCR5 rule for 3 BBAs. 
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Figure 70. The values of the resulting belief mass for 
scenario 4 and the PCR5 rule for 3 BBAs. 
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Figure 72. The values of the resulting belief mass for 
scenario 6 and the PCRS rule for 3 BBAs. 
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For the PCR5 rule for 3 BBAs, the application of the decision thresholds at the belief mass level 
ma=0,42 for scenario 4 allows us to properly evaluate the identification of the recognized object. 

For the PCR5 rule for 3 BBAs, the application of the decision thresholds at the belief mass level 
ma=0.37 for scenarios 5 and 6 allows us to properly evaluate the identification of the recognized object 
for most time moments. 


The PCR6 rule for 3 BBAs 


The simulation results of the identification information fusion using the PCR6 rule for 3 BBAs 
for the deterministic scenarios 1-6 are presented in Figures 67-72. 

For the PCR6 rule for 3 BBAs, the application of the decision thresholds at the belief mass level 
ma=0,45 for scenario 4 allows us to properly evaluate the identification of the recognized object. 

For the PCR6 rule for 3 BBAs, the application of the decision thresholds at the belief mass level 
ma=0.37 for scenario 5 and mo=0.4 for scenario 6 allows us to properly evaluate the identification of 
the recognized object for most time moments. 

Comparing Figures 64-66 with Figures 70-72 and 76-78, one can conclude that the PCR5 for 3 
BBAs and PCR6 for 3 BBAs rules provide more stable results of combined belief masses (smaller 
amplitude of changes). Due to the large dispersion of belief mass changes for scenarios 5 and 6, it is 
not possible to correctly evaluate the identification of the recognized object for all time moments. 

The presented results (Figures 61-78) allow us to conclude that the applied methods of removing 
conflicts in the information fusion allow to draw correct conclusions about the real identification of 
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Figure 73. The values of the resulting belief mass for Figure 74. The values of the resulting belief mass for 


scenario 1 and the PCR6 rule for 3 BBAs. scenario 2 and the PCR6 rule for 3 BBAs. 
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Figure 75. The values of the resulting belief mass for Figure 76. The values of the resulting belief mass for 
scenario 3 and the PCR6 rule for 3 BBAs. scenario 4 and the PCR6 rule for 3 BBAs. 
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Figure 77. The values of the resulting belief mass for 
scenario 5 and the PCR6 rule for 3 BBAs. 


Figure 78. The values of the resulting belief mass for 
scenario 6 and the PCR6 rule for 3 BBAs. 


8. Conclusions 


The proposed basic belief assignment model for ELINT-ESM sensors and radars can be used to 
build identification information fusion systems. Practical meaning have, first of all, models 
conformable to STANAG 1241. 

Due to the assumption of conflicts between the ELINT-ESM sensor declarations in this work, 
Dezert-Smarandache theory is used to determine the basic belief assignment of declarations being 
the product of the process of fusion of identification information sent by these sensors. 
Supplementing standard reports on detected signals with random identification declarations allows 
the use of methods of identification information fusion in the information fusion center. The test 
results confirm the full usefulness of conflict redistribution rules in reports from ELINT-ESM sensors 
developed as a part of Dezert - Smarandache theory, with the best results presented in the PCR5 and 
PCR 6 rule. 

The basic belief assignment model for ESM sensors and for combined primary and secondary 
(IFF) surveillance radars [15] can be applied to build models of different identification data fusion 
systems. The practical significance has first of all models compatible with STANAG 1241. It contains 
definitions corresponding to intersections of basic identification declarations. Therefore, the paper 
uses Dezert-Smarandache theory for calculation the basic belief assignment. 

The conducted research showed that the best results obtained for the PCR6 rule when reports 
from three sources (from two sensors and the fusion system database) were processed 
simultaneously. This corresponds to synchronous processing of reports and involves delayed 
processing of a report from one of the sources. The research confirmed a slight advantage of the PCR6 
rule over the PCR5 rule. This was mainly the case when the sensors sent information with a high 
degree of conflict. 
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